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Role of interactions in ^'^Rb-^'^K Bose-Fermi mixtures in a 3d optical lattice 
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We investigate the effect of interspecies interaction on a degenerate mixture of bosonic *^Rb and 
fermionic ''"K atoms in a three-dimensional optical lattice potential. Using a Feshbach resonance, 
the *^Rb-*''K interaction is tuned over a wide range. Through an analysis of the *^Rb momentum 
distribution, we find a pronounced asymmetry between strong repulsion and strong attraction. In 
the latter case, the Bose-Hubbard parameters are renormalized due to self-trapping, leading to a 
marked shift in the superfluid to Mott insulator transition with increasing Bose-Fermi interaction. 

PACS numbers: 03.75.Ss, 34.50.-s, 37.10.Jk, 71.10.Fd 
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Mixtures of quantum gases in optical lattices form 
novel quantum many-body systems, whose properties 
are governed by the interplay of quantum statistics, 
inter- and intraspecies interactions, as well as the rel- 
ative atom numbers of the components. In particular, 
degenerate Bose-Fermi mixtures, which rarely occur in 
nature and have only recently come within experimen- 
tal reach, have stimulated theoretical investigations. Al- 
ready for the simplest description of such an interacting 
mixture in terms of a single-band Bose-Fermi Hubbard 
model, a variety of quantum phases have been predicted 
including charge-density- waves, polaron-like quasi- 
particles |3|, as well as perturbed Mott-insulating (MI) 
states 0, m . The correlated (or anti-correlated) move- 
ment of composite particles can result in exotic superflu- 
ids 0. Moreover, within a limited parameter regime, 
supersolid behavior has been predicted 0, i]. Unlike 
in the pure bosonic case, realistic descriptions of the 
mixture need to consider multi-band physics finite- 
temperature [lo| . physical loss channels, and the inhomo- 
geneity present in optical lattice experiments. Recently, 
mixtures of bosonic *^Rb and fermionic ^"K atoms have 
been investigated at fixed attractive background scatter- 
ing length in two independent experiments 11 , l2| • The 



authors found a shift of the visibility loss of the °' Rb in- 
terference pattern towards shallower lattices, which has 
been interpreted as a shift of the bosonic superfluid to 
MI transition. Several mechanisms have been put for- 
ward to explain this shift, however, none could clearly be 
identified. 

Here, we report on the realization of a quantum degen- 
erate mixture of bosonic *^Rb and fermionic ^"K atoms 
with tunable interspecies interaction in a 3d optical lat- 
tice. We investigate the coherence properties of the 
bosonic component in the presence of fermions over a 
wide range of interspecies interactions, from strongly at- 
tractive to strongly repulsive, for varying admixture of 
fermions. At large fermion fillings and attractive inter- 
species interactions, we find good agreement with a re- 
cent model of self-trapping [§] . This leads to a significant 
renormalization of the Bose-Hubbard parameters, result- 



ing in a shift of the superfluid to MI, in good agreement 
with our experimental results. 

Starting point for our experiments is a degenerate mix- 
ture of 4 X 10^ ^''Rb and up to 3 x 10^ '^''K atoms in their 
respective hyperfine ground states |_F = \,mF = +1) 
and \F = |,to_f — — §)• The mixture is held in a hori- 
zontally crossed pancake-shaped dipole trap, operated at 
a wavelength of 1030 nm. The ellipticity of the dipole 
beams makes the trap vertically tight enough to ensure 
interspecies overlap in the presence of gravitational sag. 

A homogeneous magnetic field is applied in order 




-1.0 



-0.5 



I 



-200 , , 

-400 abfiao) 



• • • 



FIG. 1: Dependence of bosonic interference pattern on in- 
terspecies interaction, (a) Visibility of the *^Rb interfer- 
ence pattern versus lattice depth Vo and interaction strength, 
parametrized by the interspecies scattering length aj,/, for 
high Nj (interpolation from 1004 individual measurements). 
(b)-(e) Exemplary TOF images of *^Rb atoms released from 
a 9-Er deep optical lattice. The values for the interspecies 
scattering length are (b) abj = — 160ao, (c) abj = — 80ao, (d) 
ahf — -1-5 ao, and (e) aj,/ = -|-235ao. 
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to address the 2.9 G wide interspecies Feshbach reso- 
nance (FBR) at 546.9 G We can continuously tune 
the interspecies scattering length ai,f between — ITOoq 
and +800 ao below the FBR, and between — SOOoq and 
—200 ao above the FBR, with an accuracy of ±10 ag close 
to the resonance, where ao is the Bohr radius. After 50 ms 
settling time for the magnetic field, a 3d optical lattice 
is adiabatically ramped up within 100 ms, reaching a fi- 
nal depth in the range of Vo = 2 to 17 Er. The lattice 
wavelength A — 755 nm is blue-detuned with respect to 
the D2 transition of both ^^Rb and "^^K and chosen such 
that the lattice depths, as measured in units of the re- 
spective recoil energies = /{2mi,j X^), are equal 
for both species. This ensures maximum overlap of the 
Wannier functions, while the ratio of the tunneling rate 
of bosons and fermions is given by ruf /mi,. The blue lat- 
tice also causes an underlying harmonic anti-confinement 
which is (over-)compensated by the dipole trap (Hj . 

After a hold time of 100 ms, which is long compared to 
the tunneling time for any lattice depth used, all traps as 
well as the magnetic field are instantaneously switched 
off, and the atom clouds are allowed to expand during 
18 ms time of flight (TOF). The resulting interference 
pattern of the ®^Rb atoms is recorded using standard ab- 
sorption imaging [see Fig. [T]^b)-(e)]. From these images, 
we extract the contrast in terms of the visibility [l4l|. 
Furthermore, we determine the condensate fraction from 
a fit to the interference pattern. 

By varying both the lattice depth and the magnetic 
field, we probe the phases of the mixture. We have an- 
alyzed scenarios with low {Nf w 0.25iVf,), intermediate 
(Nf R:! 0.5iVb) and high {Nf « 0.75Nb) atom num- 
bers. The visibility diagram of the mixture with high Nf 
is displayed in Fig. [TJ Profiles of the same data are shown 
in Fig. [21 illustrating the most prominent features. 

For shallow lattices of less than 3Er, we find a high 
visibility, which is almost independent of the interaction 
strength. Also for deeper lattices, at atf ^ 0, we recover 
a monotonic decay of the visibility versus lattice depth 
compatible with the superfluid to MI transition in a pure 
^^Rb sample [l3|. This demonstrates that the fermion 
cloud becomes fully transparent for the bosons. When 
the absolute value of the scattering length is increased, 
the visibility decay is shifted towards shallower lattices. 
For small absolute values oiabf, the decrease of visibility 
is symmetrical around a^f ~ 0. The symmetry persists 
within the interval \Ubf\ ^ Ubb^ where the interboson 
repulsion is strong enough to effectively hinder bosonic 
double occupation of sites. Departing from this symme- 
try peak, we see a further decay of visibility, which is 
significantly stronger for intermediate attractive than for 
comparable repulsive scattering lengths. This suggests a 
fundamental difference in the underlying mechanisms on 
either side. The behavior in the vicinity of the attractive 
background scattering length as we observe it is similar 
to earlier measurements ll|, |12| . Towards very strong at- 



traction, we find a significant loss of ®^Rb atoms, accom- 
panied by an increase of visibility. For strongly repulsive 
interactions, the visibility remains almost constant on 
a high level. This observation is consistent with phase 
separation in the lattice between bosons and fermions, 
which might be promoted by the slightly different trap 
shapes they experience, or the formation of a strongly 
anti-correlated mixed phase. Presently, we can not dis- 
tinguish between both scenarios. For abf > 400 oq, we 
observe significant atom losses, however, these are es- 
sentially independent of the hold time, suggesting their 
occurance early during the lattice loading, when the mix- 
ture is held just below the FBR. 

Based on the observations outlined above, we quali- 
tatively divide the diagram of Fig. [1] into five distinct 
regimes, i.e. lossless coexistence in very shallow lattices 
for all scattering lengths, and regimes characterized by 
very strong interspecies attraction, intermediate attrac- 
tion, weak interaction (of either sign), and very strong re- 
pulsion for deeper lattices. This classification also holds 
for low and intermediate values of Nf , although the sym- 
metric feature around abf = is most pronounced for the 
highest ''^"K numbers. 

We have verified that the visibility decrease is re- 
versible for interspecies scattering lengths in the range 
— 200ao < abf < 400 ao, when abf is tuned back to zero. 
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FIG. 2: Profiles of visibility versus interspecies scattering 
length at lattice depths of (a) 3Er, (b) 5Er, (c) 9Er and (d) 
12 Er for high Nf. At all lattice depths beyond 3 Er, the visi- 
bility shows a maximum at a position consistent with ai,f = 
as indicated by the arrow within the calibration uncertainty. 
Blue diamonds (red circles) indicate points measured above 
(below) resonance, as depicted in the inset of (a). The shaded 
region represents the visibility measured for a pure *^Rb cloud 
using the same experimental parameters. 
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FIG. 3: Onset of the loss-dominated regime, (a) Experi- 
mentally determined visibility minimum and (b) simultane- 
ous *'^Rb loss feature at Vo = 9i?r for intermediate Nf. (c) 
Critical ai,/ determined from the visibility data, versus lat- 
tice depth for high, intermediate and low Nf (darker colors 
for higher Nf). The shaded area gives an estimate of the scat- 
tering length at which the bosonic three-body loss timescale 
becomes shorter than the experimental hold time, based on 
our variational model. The inset shows a typical TOF im- 
age of *^Rb atoms in the loss-dominated regime, for low Nf . 
The cross indicates the associated experimental parameters 
[Vo = 11 Er, Ubf = -480 ao). 



provided that the system is given enough time to reequi- 
librate. This suggests that the reduction of coherence is 
caused by a reversible redistribution of the atomic densi- 
ties during the lattice loading. For attractive scattering 
lengths beyond a critical value, we observe significant 
losses of ^^Rb atoms, along with an increase in visibil- 
ity. Despite this increase, no condensate is discernible 
in the corresponding TOF images. These findings are 
shown exemplarily in Fig. EJa), (b). We fit two straight 
lines to the relevant part of the visibility profiles, and 
use their intersection point as an estimate for the critical 
value. Fig. [2{c) shows that the resulting critical scatter- 
ing lengths are independent of Nf. This suggests that the 
visibility increase at strong attractions is governed by the 
onsite physics. Recently, the crucial role of the modifica- 
tion of the Wannier orbitals due to the mean-field onsite 
interaction in understanding the attractively interacting 
mixture has been pointed out 0. A strong interspecies 
attraction will not only enforce high boson occupation 
numbers, but will also at given occupation yield increased 
onsite densities due to admixture of higher bands, which 
in turn could lead to a significant increase in three-body 
loss rates. We estimate this onsite density increase from 
a variational harmonic oscillator model, minimizing the 
total onsite energy with respect to the width o-f, and a/ 
of the bosonic and fermionic Gaussian density profiles 
for given lattice depth and occupation numbers. For 
all practically relevant boson numbers, we find a critical 
scattering length at which the onsite density collapses. 
The strong density increase in the vicinity of this col- 
lapse gives rise to an enhanced bosonic three-body loss 



rate due to the N3 oc cr^ scaling. For practical pur- 
poses, a site can be considered lost when N^t ^ 1, where 
T is the experimental hold time. Assuming a loss coef- 
ficient for collisions of three ®^Rb atoms in the range 
K3=4...7 X 10~^° cm*^ /s [11], we find that lattice sites 
with four or more bosons will be lost on time scales faster 
than our hold time even in the absence of fermions. On 
the other hand, the observed ratio of lost ®''Rb to 
atoms between 3 and 4 suggests that three-body losses 
of two bosons together with one fermion do not play a 
significant role in our system. This observation is con- 
sistent with the measurements presented in [llj. There- 
fore, we focus our attention on sites occupied by three 
bosons and one fermion. The critical scattering length 
for bosonic three-body losses derived from our variational 
model shows good agreement with the experimentally ob- 
served onset of the loss-dominated regime, as can be seen 
in Fig.[3l^c). We therefore conclude that for strong attrac- 
tion, the dominant process is the loss of highly occupied 
sites enhanced by the interaction-induced density accu- 
mulation, commonly referred to as self-trapping . 
The associated increase of the visibility could be thought 
of as the removal of strongly localized atoms on highly 
occupied sites from the system, and thus a reduction of 
the bosonic filling in the system. 

In order to investigate the persistence of bosonic super- 
fluidity in the mixture, we extract the condensate frac- 
tion from the recorded TOF images. This is done by fit- 
ting anisotropic two-dimensional bimodal Gaussians to 
all first-order diffraction features as well as the central 
peak after substraction of a broad Gaussian background 
from the image. The condensate fraction, which we de- 
fine as the sum of the atom numbers in the narrow fea- 
tures at quasi-momentum zero and the eight first-order 
peaks, divided by the total atom number 2M ' found to 
decay monotonically towards zero with increasing lattice 
depth for all values of abf at all fermion numbers. 

The point in lattice depth at which the condensate 
fraction vanishes, depends strongly on the interaction 
strength (Fig. 2]). To quantify this point, we fit a lin- 
ear decrease model to our data, thus obtaining the kink 
position for each scattering length (inset of Fig. 2]). For 
medium to high Nf, the transition points lie almost on 
top of each other, both displaying a strong shift of up 
to 10 Er towards shallower lattices for attractive interac- 
tions. In this regime, we expect the bosons to move on 
top of an essentially homogeneous fermion filling. The 
fermions then constitute an effective attractive potential 
for the bosons, which adds to the depth of the lattice 
potential. The shift can thus be understood in terms 
of an effective deepening of the lattice sites induced by 
the Bose-Fermi attraction, which modifies both the bo- 
son tunneling and interbosonic repulsion. As soon as 
the interspecies attraction overcomes the bosonic repul- 
sion, the formation of a bosonic MI shell of filling two 
is energetically favored, resulting in even deeper effective 
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FIG. 4: Shift of MI transition point with interspecies interac- 
tion. The diamonds and circles represent experimental runs 
at medium and high fermion numbers, respectively. The error 
bars incorporate the fit uncertainty and the weighted average 
over typically two experimental runs. The lines indicate the 
theoretical MI transition point as estimated by the variational 
model for filling one (solid gray line), as well as exact diag- 
onalization predictions for MI fillings of two (solid red) and 
one (solid green). For atf > 0, the solid lines correspond to a 
phase separation scenario, while the respective mixed phase 
predictions are depicted by dashed lines. The inset exemplar- 
ily shows the behavior of the condensate fraction at interme- 
diate Nf, for ttbf = — 295ao (black points) and abf — 235 ao 
(green points), respectively. 



potentials, thereby further enhancing the modification of 
the Hubbard parameters. The variational model (gray 
line in Fig. |4]) clearly underestimates the self-trapping 
and thus yields smaller shifts than experimentally ob- 
served. A single-site exact diagonalization approach for 
the effective Hubbard parameters including the full or- 
bital degrees of freedom ^] (solid red and green lines 
in Fig. [3]) shows good agreement with the experimental 
data This supports that the shift of the transi- 

tion point in the attractive regime can be explained by 
self-trapping. For repulsive interactions, a shift towards 
deeper lattices of comparable magnitude is predicted by 
the same model, assuming a homogeneous distribution 
of bosons and fermions throughout the lattice. The ab- 
sence of a significant shift in this direction indicates that 
bosons and fermions do not occupy the same lattice sites 
in that regime. For low Nf, where the assumption of a 
homogeneous fermion filling is not valid, the decrease of 
the condensate fraction shows a more complex behavior 
which can not be captured by a linear slope, and the 
discussion of which is beyond the scope of our model. 

In conclusion, we have explored the quantum phases of 
a mixture of bosonic ®^Rb and fermionic '^''K in a 3d opti- 
cal lattice. By tuning the interspecies interaction and lat- 
tice depth, we have observed transparency of the fermions 
at ahf — 0, and a symmetric decay of visibility in a small 
window around that point. For stronger interactions, we 
find a remarkable asymmetry between the attractive and 



repulsive side. The condensate fraction shows a behavior 
which we interpret as a shift of the bosonic MI transition 
towards shallower lattice depths with increasing inter- 
species attraction. We attribute this shift to a renormal- 
ization of the parameters of the effective Bose-Hubbard 
model by self-trapping, and find good quantitative agree- 
ment with an exact diagonalization calculation. Our re- 
sults suggest that Feshbach-assisted lattice loading can 
be used to tailor the resulting lattice filling. Finally, the 
observed transparency is a promising starting point for 
experiments on tunable impurities 17, 18| . 
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